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SUMMARY 

An effective procedure for NASTRAN has been developed that permits any number of substructures 
bf any size to be synthesized for the purpose of developing normal modes of vibration of the complete 
structural system. The technique is extended to permit modal transient analysis of the subdivided 
bystem. This latter procedure permits the use of NASTRAN’s ability to include nonlinear forces in the 
problem. The five-phase process is accomplished using standard ttASTRAN rigid formats with 
problem-independent alter packages and DMAP sequences. 

I 

1 INTRODUCTION 

(The advent of finite-element techniques has led to development of computer programs that perform 
efficient analyses of large-order structural systems of arbitrary geometry. Yet, for very large problems, 
it becomes impractical to attempt to develop and analyze a mathematical model of the entire 
structure. Questions of adequate data management and the cost of the analyses make such efforts 
difficult. Consequently, techniques have been developed to permit the structural system to be divided 
into smaller units, or substructures, which may then be modeled separately. 

Substructure models, represented by some truncated set of degrees of freedom, are eventually 
brought together by the application of compatibility requirements at the interfaces. The resultant 
system representation is analyzed and the behavior of each substructure is determined by a 
back-substitution process. 

The mode synthesis technique has been developed to permit dynamic analysis of such subdivided 
structures. A truncated set of vibration modes of a component is used as generalized coordinates to 
approximate dynamic characteristics. The several substructures are coupled by ( 1 ) compatibility 
relationships at the interfaces, and (2) modal transformations relating physical degrees of freedom 
with the generalized component modal coordinates. Solution of the resultant equations yields system 
eigenvalues and eigenvectors. 

f The use of vibration modes alone, however, usually fails to describe adequately the behavior of 
the system, particularly in the regions of the interior boundaries (interfaces). Consequently, the 
[substructure vibration modes are augmented by other displacement functions. Numerous alternatives 
will be found in the literature. Hurty (Reference 1) and Craig and Bampton (Reference 2) employ 
“constraint” modes, defined by unit displacements of interface degrees of freedom on a fixed 
interface boundary. Bamford, et al (Reference 3) and Rubin (Referenced) use “attachment” modes, 
defined by unit forces at each interface degree of freedom. Rubin has shown that this concept may be 
used for free interface boundaries, as well as fixed, to produce inertia relief deformations in 
unrestrained substructures. It is this latter arrangement that is employed in the development of mode 
synthesis via NASTRAN that follows. The technique can be easily modified to use constraint modes 
and fixed interface boundaries. 


Synthesizing substructures for dynamic analysis via NASTRAN has been discussed previously. 
The NASTRAN Theoretical Manual (Reference 5) discusses the use of scalar elements and MPC sets to: 
represent substructures. These procedures were reviewed and amplified by Courtney (Reference 6). 
The current approach exceeds this earlier work by explicitly including attachment modes and by 
providing automated communication of data between component and system levels. The ability to, 
develop the structure characteristics and acquire results in spatial coordinates enhances the procedure,; 
particularly when considered from a production environment viewpoint. 

The synthesis may be performed using any number of independently formulated substructures of, 
any size. No restrictions are placed on the identification schemes used in component model 
development, since they are not carried over into the synthesis. The technique uses standarized alter 
packages and DMAP sequences that are virtually independent of the problem characteristics! thus 
the often present necessity for user specification of many data handling directives is avoided . Control 
for matrix partitioning and merging are in the form of user-supplied partition vectors. I 

An extension of the mode synthesis procedure includes the modal transient response analysis; 
for the synthesized system, wherein the required modal data is recovered on a substructure- j 
by-substructure basis. Since the complete system may be quite large, the modal matrix is i 
truncated to include only those degrees of freedom subjected to external loads or required to definej 
the presence of nonlinear forces. Such nonlinear forces are included only after the mode synthesis has 
been completed. The responses of the system, including internal forces and stresses, if desired, are] 
recovered by back-substitution into each substructure. , 

The complete analysis process is divided into five sequential phases, with communication , 
between phases accomplished via User Tapes; e.g., INPT, and Problem Tapes (NPTP, OPTP). 

SUBSTRUCTURE MODAL ANALYSIS 

Phase I of the mode synthesis/transient response analysis using NASTRAN is the development of 
substructure, or component, modal data. This information is assembled in Phase II for system modal 
analysis. There are few restrictions and very little additional demands required for this phase over that 
associated with a standard NASTRAN analysis. It is necessary to give special consideration to those 
degrees of freedom in the component model that are involved in connection with other portions of the 
structural system. Additional information is required to direct the preparation of the component 
modal data for later synthesis. 

The full range of NASTRAN’s resources may be employed to define a substructure model. Since 
substructures are independently defined, no special restrictions are placed on identification or 
numbering schemes. 

Components may have any form of boundary condition at their external boundaries; i.e., those 
boundary grid points not connected to other components. Connection grid points must be 
unrestrained in those coordinates that are used to join the components. 

Since the interface degrees of freedom are, by definition, unrestrained, the collectivity of normal 
vibration modes cannot correctly define internal loads at the interface of the component. To correct 
this deficiency, the dynamic modes are augmented by a selected set of static deformations, which, 
when properly weighted and summed, will provide the internal loads at the interface degrees of 
freedom. The static modes are generated by applying a unit force at each interface degree of freedom 
in turn and letting the component deform. If rigid body motion is prevented by external constraints, a 
straightforward static analysis results. If, however, a rigid body component of motion does result from 
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the unit load at the interface, an inertia relief analysis is required to produce internal loads. For this 
jreason, Phase 1 of the mode synthesis procedure is performed using Rigid Format No. 2, Static 
[Analysis with Inertia Relief, as modified by the alter package of Figure 1. 


4 PHASE ONE ALTER PACKAGE 

ALTER 

so* so 

PURGE 

KRR»KLR*QR*DM/REACT 4 

CHKPNT 

KHR»KLR»0R»DM 4 

alter 

76 

EOUIV 

KAA»KLL/REACT 4 

CHKPNT 

KLL 4 

ICGND 

LBL6* REACT $ 

ALTER 

7B 

LABEL 

LBL6 4 

alter 

80 

COND 

LBL7»REACT 4 

alter 

84 

! LABEL 

LbL7 4 

1 ALTER 

88 

| EuUIV 

pl»pli/rlact/po*poi/react 

! COND 

LbLIO* REACT 4 

! ALTER 

89 

label 

LBLlu 4 

1 ALTER 

110,* 114 


4 VIBRATION modes 

DPD UYNAMICS»GPL»SIL#USET/GPLD'SU.D»USETD» » » » * » »EED* EQcYN/ 

V,N,LUSET/V,N,LUSETD/V,N,^ 1 /V,N,Z2/V»N.25/V,N,Z4/ 
V*N»ZS/V*N*Z6/V»N*Z7/C*N*G/V*N»Z9 4 

l SAVE Z9 * 

, PARAM //C,N, ADD/V,N,VI6CASE/V, Y»NCON=0/C,N,1 4 MODES SUBCASE COUNT 

Read i\aa»maa*mr*dm*eed*uset,casecc/lama»phia*mi,oeigs/ 

C *N» MODES/V * N* NEIGS/V * N* V IBCASE 4 
SAVE NEIGS * 

CHKPNT LAMA,PHIA,MI,OEIGS * 

OFP LAMA,0EIGS* * , » //V » N* CARDNO $ 

SAVE CARDNO 4 

COND FINIS* NEIGS 4 

* RECOVER FULL 'MODE VECTOR 

SuRl USt.T » * PHI A **»GO*GM*»KFS»» /PHIG » * GGX/C *N»1/C»N*REIG 4 

CHKPNT PHIG » GGX $ 

SDR2 CASECC*CSTMrMPT*DIT *EQEXIN,SIL* * *BGPDT * LAMA »QGX* PHI G*EST * / 

* 0uG2 * OPHIG* 0ES2 * 0EF2 »PPHIG/C »N* REIG 5 
CHKPNT PPHIG 4 

OFP UPHlG* OQG2 * 0EF2 » 0ES2* * //V * N» CARDNO 4 

SAVE CARDNO $ 

4 DYNAMICS PORTION OF GENERALIZED MATRICES 
i GKAM» *PHIA*MI ,LAMAr * » * »CASECC/MHH»BHH*KHH»PHH/V»N»Z9/ 

C,N,iOO/C,N,0./C,N*0./C*N,-.1/C»N,-i/C»N,-:l/V»N,NOC/V,N,FMODE 4 
4 ATTACHMENT MODE PARTITION OF GENERALIZED MASS AND STIFFNESS MATRICES 
MRYAD UGV*MG6*/AM/C*N»1 4 

MRYAD AM,UGV*/MGEN/C,N,0 4 GENERALIZED MASS (ATTACHMENT) 

MpYAD UGV * KGG* / AK/C * N* 1 4 

MPYAD AK»UGV»/KGEN/C*N»0 4 GENERALIZED STIFFNESS (ATTACHMENT) 

Figure 1 . Phase I alter package. 
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CHKPNT MGEN , KGEN % 

% OKF-OIAGONaL partitions 

MPYAO AM.PHIG»/ML/C*N»0 5 L.L. PARTITION 

MPYAD AK,PHIG./KL/C,N»0 * L.L. PARTITION 

TRNSP ML/MO S U.R. PARTITION (MASS) 


(MASS) 

(STIFFNESS) 


TRNSP KL/KU * U.R. PARTITION (STIFFNESS) 

S GENERALIZED MATRICES AND CONNECTION PARTITION OF MODES 
MERGE MHH,ML»MU,MGEN t STATlC»/MASS/C,N»-l/C»Nr2/C»N,6 * 

MfcRGE KHH*KL*KU*KGEN*STATIC*/STIFF/C*N*-1/C»N»2/C*N*6 * 

MERGE PHIG* *UGV* »STATIC»/M0DES/C#N»1/C»N,2/C#N»2 $ 

PARTN MODES* *COOF/*INTER» */C*N,l/C»N»2/C*NrO/C*N*2 * 

MERGE GGX* *QG* *STATIC*/QQG/C»N»1/C»N*2/C*N*2 S 

ChKPNT MODES »GQG * 

PARAM //C*N*NOP/V*Y*PRINT=-l * 

COND list*print * 

MATPRN MASS, STIFF* INTER* *//* 

LABEL LIST 5 

* WRITE MATRICES ONTO USER TAPE -INPT- 

PARAM //c»N*SUb/V*N*KOUNT/V*Y*COMP=0/C»N»1 s 

PARAM //C,N*MP r /V,N,C0UNT/V > N»K0UNT/C»N,3 * 

COND FINIS* COUNT $ EXIT IF SEQUENCE NUMBER IS LESS THAN U 
PARAM //C * N * SUB/V » N * IOTEST/ V , N , COUNT/C * N, 1 $ 

COND IOi , lOTEST S JUMP IF FIRST COMPONENT 

INPUTT1 /»***/C*N*-l * REWIND -INPT- 

JOMP 102 * 

LABEL 101 $ 

output i » »»*,//c»n,-i » 
label ioz * 

OUTPUT 1 MASS* ST IFF * INTER* *//V *N*COUNT S 
ENDALTER 


Figure 1. Phase I alter package ( concluded). 

The DMAP alterations supplied for Phase I add the eigenvalue extraction module and other 
modules that merge the results of the static and dynamic processes to form the final component mode 
matrix. 

Data required by Phase II is output to a User Tape. These data include the generalized mass and 
generalized stiffness matrices associated with the merged mode vector matrix and the merged mode 
vector matrix truncated to include only the interface degrees of freedom. The truncation is directed 
by a partitioning vector supplied by the user. The Phase I analysis should be checkpointed for later use 
in recovery activity in Phases III and V. 

Additional data required beyond the model description include special control parameters 
entered via bulk data card PARAM, and partition vectors, entered via bulk data card DMI. The specia! 
parameters required for Phase I include: 

1. COMP — The sequence number of the component. Used to position the User Tape before writing 
component data. 

2. NCON — The number of interface degrees of freedom for the component. The number of stati* 
load subcases/static modes must agree with NCON. 


3. PRINT — A positive value causes printing of the matrices written on the User Tape. 


plie partition vectors needed are: 

jl . STATIC - Used to merge static and dynamic modes into one mode vector matrix. M = sum of all 
modes computed (static and dynamic, including rigid body). An element value of 1 locates astatic 
mode in the merged matrix. There must be exactly the same number of positive elements in the 
| vector as there are static modes. 

2. CDOF - Used to identify the interface degrees of freedom. M = g size. Vector contains 1 for 
elements corresponding to the internal sequence number of all interface degrees of freedom. 

j The Case Control deck must include (n + 1 ) subcases, where n is the number of interface degrees 
of freedom for the component (and, consequently, the number of static load cases). The 
first n subcases are static cases, with a LOAD reference in each. The last subcase is the modal 
extraction subcase, with a METHOD card. Figure 2 illustrates a typical Case Control Deck. 

REQUEST(NPTP,HD»S) ring in 
j REQUEST (INPT# HD# S) RING IN 

ATTACH! Nib »NASTRAN1S#ID=TONE90970) 
j REMIND (NPTP» INPT) 

MAP (OFF) 

N15. ATTACH 
RETURN (NPTP, I NPT) 

EXIT. 

DMPU40000) 

RETURN(NPTP,INPT> 

(7-8-9) 

ID PHASE# ONE 
APP DISP 

SOL 2# 1 

ChKPNT YES 

time 2 

s (PHASE one ALTER PACKAGE GOES HERE) 

CEND 

T1TLE=DEM0NSTKATI0N OF PHASE ONE DECK SET-UP 
SPC = 173 
MPC = 4177 

output 

SET 37 = 20I THRU 287 
VECTOR = ALL 
SPCFORCES = 37 

* THIS EXAMPLE ASSUMES FOUR CONNECTION 

* DEGREES OF FREEDOM. 

subcase 100 

LOAD =57 
SUBCASE 200 
LOAD = 371 

Subcase 300 

LOAD = bb 

subcase 4uu 

LOAD = 331 
SUBCASE 500 
METHOD = BbO 

begin bulk 

* (SULK DATA DECK GOES HERE) 

enddata 

( 6— 7-8-9) 

Figure 2. Typical deck setup for Phase I. 
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MODE SYNTHESIS 


Phase II of the mode synthesis analysis involves generation of system dynamic characteristics throu, 
synthesis of component information. This is performed using Rigid Format 3, Normal Mode Analys 
modified by the alter package of Figure 3. Direction is required from the user regarding compatibl y 
relations among the various components, and selection of component generalized coordinates. . 
attention to certain aspects of internal data storage arrangements is required; in most cases, howev, 
this is a trivial consideration. Since the activity involves operations on generalized, rather th 
physical, coordinates, the intuition and understanding of the user become much more necessary th 
in a normal study. Further, component coupling is accomplished by imposing transforms (develop 
from physical coordinate relationships) on the generalized coordinates. This requires the use of sor 
mathematical sleight-of-hand that can lead to confusion unless the user is alert. 


$ PHASE TWO ALTER PACKAGE 

^Itore' component PARTITIONING vectors on temporary disk file 
0UTPUT2 P0S1 * LOCI * TRUNK 1 * *//C*N*“l/V» Y *HOLD = H * 

0UTPUT2 P0S2*L0C2*TRUNK2* »//C*N»0/V*Y»HOLD * 

$ INSERT ADDITIONAL -0UTPUT2- STATEMENTS 
* UNTIL ALL COMPONENTS ARE TREATED 
INPUTT2 /*»»»/C*N*-l/V*Y» HOLD » 

PARAM //C*N*NOP/V*N,TRUE=-l S 
% PARAMETER K IS COMPONENT LOOP COUNTER 


PARAM 
PURGE 
INPUTT1 
JUMP 
INPUTT2 
$ LOOP 
$ merging 

LABEL 
INPUTT 2 
INPUTT1 

partn 
partn 

PaRTN 

merge* 
add 

EUUIV 

merge* 
add 
eguiv 
merge* 
add 

EuUIV 
PARAM 
COND 
REPT 
LABEL 
CHKPNT 
PARAM 


//C*N*SUB/V*N*K/V* Y*NC0MP=1/C,N»1 S 
SUMK*SUMM,SUMC/TRUE * 

/* * * */C*N*-l $ 

TOPI * 

/SUMK*SUMM*SUMC* */C*N*0/V»Y»XTRA=13 * 

OVER ALL SUBSTRUCTURES, TRUNCATING COMPONENT DATA AND 
TO FORM UNCOUPLED SYSTEM MATRICES 

/cp1rp*trunk**/c*n*o/v*y,hold * partition vectors 

/MASS* STIFF, INTER* ,/C *N»0 * COMPONENT DATA 
MASS* TRUNK* /MM* * ,/C*N*— 1/C»N»2/C*N,6 * 

STIFF *TRUNK, /KK *, ,/C, N,-l/C,N, 2/C *N,t> * 

INTER * TRUNK */CC* * * /C * N, 1/C »N* 2/C * N* 2 * 


$ 


, *CC*CP*RP/CCX/C*N*1/C»N*2/C»N»2 
SUMC*CCX/CCXX S 
CCxX * SUMC/TRUE % 

* * * MM * CP* /MMX/ C *N»-1/C»N,2/C»N*6 
SUMM, MMX/MMXX * 

MMXX * SUMM/TRUE * „ „ 

* * , KK * CP * /KKX/C *N,-1/C»N»2/C»N»6 % 

SUMK * KKX/KKXX * 

KKXX,SUMK/TRUE % 

//C*N»SUB/V*N*K/V *N*K/C*N*1 SDECREMENT COUNTER 
END1*K * 

TOPI *20 $ 

ENU1 * 

SUMC, SUMM, SUMK * 

//C,N*NOP/V*Y*PRINT=-l * 


Figure 3. Phase II alter package. 
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COND PRNTl » print 4 

MaTPRN SUMC t SUMM» SUMK » > //4 PRINT SYSTEM UNCOUPLED MATRICES 
LABEL PKNT1 * 

i PARAM //C»N»NOP/V»N.NSKIP=0 $ 

GP4 CASECC » GE0M4 ,EQEXIN* SIL » GPDT /RGG» # USET » / 

I V »N» LUSE T/V » N» MPCFl/ V » N» MPCF2/V»N» S INGLE/ V » N» OMI T/ 

V »N» REACT/ V» N# N5K IP/V » N» REPEAT /V t N» NOSET/V»N»NOA/ V » N » NOL 4 

save mpcfi »mpcf2* single* omit # react ,nsk ip » repeat » noset »noa,nol s 

CHKPNT KGG, USET $ 

* compatibility relations yield reduction transform 
vec uset/vf/c»n»g/c»n#comp/c»n»f s 

PARTN RGG»VF#/RG# » */C»N»l/C»N»2/C>N»2 4 

I MPYAD RG»SuMC»/A/C*N»0 4 COMPATIBILITY IN GENERALIZED COORDS. 
VEC USET/VO/C ,N» F/C »N, COMP/C , N,0 4 

PARTN A»VO»/AI » » AD»/C »N» 1/C »N*2/C »N» 2/C»N* 0/C »N» 1 4 

i SOLVE AD»AI/B/C#N»0/C»N»-1 4 

| COND PKNT2* PRINT 4 

, MATPRN A#b» * * //4 

! LABEL PKNT2 * 

I SMP2 USET»B»SUMM/MAA 4 REDUCED SYSTEM MASS MATRIX 

[ SMP2 USET ,b»SUMK/KAA 4 REDUCED SYSTEM STIFFNESS MATRIX 

COND PRNTi* PRINT 4 

j MATPRN MAA,KAA, , ,//$ 

LABEL PRNTS 4 

l 4 CONTINUE WITH STANDARD ANALYSIS METHODS USING -MAA- AND -K AA“ 

; alter 8o*au 

FbS LLL,ULL,KLR/DM/C,N, 1 /C»N,-i/C,N,2/C»N,2 4 

alter s»s# 112 

! MPYAD B»PHIA»/uEPEND/C»N»0 4 RECOVER DEPENDENT COORDS 

i merge phia, depend* »* # vo/full/c»n»i/c»n,2/c»n»2 s 

0UTPUT2 LAMA, MI, FULL, ,//C,N,0/V,Y»HDLD 4 

4 PRINT ON DEMAND 

PARAM //C»N#NOP/V»Y,VECTOR=-l 4 

i COND PRNT4» VECTOR 4 

MATPRN FULL , * # > //4 

LABEL PKNT4 4 

ENDALTER 

Figure 3. Phase II alter package (concluded). 


All component data required for mode synthesis is, at the beginning of Phase II, contained on a 
Jser Tape in a sequential manner. In most cases, the number of component generalized coordinates 
nodes) to be used in the synthesis task will be less than the number provided through the individual 
omponent analyses of Phase I. The user identifies the modes to be used by supplying partitioning 
lectors for each component. 

The key to the mode synthesis task is the compatibility relationship among the interface 
oordinates, expressed in the following equation: 


s x (l) + s 

1 c 2 c 


+ s.x u; = o 


/here X are the interface coordinates for component i and S. is a transformation employed to 
C 1 
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For simplicity, the compatibility equation (Eq. 3) is rewritten as 

[A] {q} =0 (4) 

Motion at the interfaces is redundantly defined in the above equation. A set of dependent coordinates 
must be identified such that they are equal in number to the excess interface coordinates, which, in 
turn, are equal in number to the compatibility equations. The choice is essentially arbitrary, excepl 
that for a singular system an appropriate set of rigid-body modes must be retained as independenl 
coordinates for reasons explained below. This set must be adequate to supply rigid-body motion in al 
directions that are singular in the complete structural system. The compatibility equation i: 
partitioned into dependent and independent coordinates: 

> 

q i 

• • • ► = 0 (5 

q D 
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( 6 ) 


rhe relationship between dependent and independent coordinates is thus 

i 

1 V - -‘V 1 'V <v 

I 

It is apparent from Eq. 6 that the dependent coordinates must be chosen so that [Ap] is 
nonsingular, and hence invertible. The full set of generalized coordinates can now be defined in terms 
of the independent coordinates. 

{q i } 

( 7 ) 



This equation supplies the necessary coupling of component generalized coordinates required to define 
the behavior of the full structural system. Letting 



( 8 ) 


the uncoupled component generalized mass and stiffness matrices are transformed into coupled system 
matrices. If the assembled uncoupled matrices are [ HI. ] and [k] where 



3^9 



then application of the transformation results in 

[M] = [/3] T [ ml [0] j 

[K] = [0] T [k] m 

where [M] and [K] are coupled system mass and stiffness matrices in terms of the independent set of, 
generalized coordinates. These are used in the eigenvalue extraction to acquire system model data. The 
mode vectors obtained are the system vibration modes, expressed in terms of the independent! 
component generalized coordinates. 

The system mode vectors are expanded to include the dependent generalized coordinates and 
stored, along with the eigenvalue table and system generalized mass matrix, on another User Tape. 
These data are accessed in Phase III for purposes of recovering mode data in component physical) 
coordinates. | 

Each component generalized coordinate is represented in the NASTRAN analysis by a scalar] 
point. In addition, the interface degrees of freedom are mapped on a one-to-one basis on another set of | 
scalar points. It is advantageous to employ a numbering sequence that can be used to distinguish easilyj 
between those scalar points representing component generalized coordinates and those representing 
interface degrees of freedom. The numbering scheme for the scalar points is arbitrary; however,! 
partitioning vectors described below must be compatible with the internal sequence of points-, thus, it ] 
is desirable that some logical numbering scheme be employed to facilitate bookkeeping. Both! 
component mode and component interface degree of freedom internal sequences cannot be altered 
from the Phase I condition. Regardless of what scalar point numbering scheme is employed, 
compatibility with these sequences must be maintained. Note that, if necessary, SEQGP cards can be 
used to resequence scalar points. 

The definition of interface compatibility relationships • (Eq. 2) is introduced by means of 
multi-point constraint equations. It is the responsibility of the user to verify the coordinate 
relationships defined by the MPC equations, including any necessary coordinate transformations. The 
MPC cards supply these relationships in terms of the scalar point equivalents of the interface degrees 
of freedom. The relationships defined must be sufficient to eliminate all redundancy in describing 
interface motion. To illustrate, if there are n components interfacing at a grid point, then there 
must be (n- 1) MPC equations supplied for every direction in which motion can occur at that grid 
point. Within that requirement, the choice of dependent degrees of freedom is arbitrary, since the 
actual synthesis activity does not consider the distinction. The scalar points used as independent 
coordinates in the MPC equations are placed in the S-set by SPC cards. This is done solely for the 
purpose of internal bookkeeping and does not impose any real constraint on the actual structure. The 
transform to generalized coordinates (Eq. 3) is performed automatically. 

The division of component generalized coordinates into independent and dependent sets, as 
required by Eq. 5, is done by specifying the dependent set on OMIT cards. The number of coordinates 
omitted must be identical to the number of interface dependency equations defined by MPCs. The 
specification is in terms of the scalar points representing the generalized coordinates. Based on this set 
selection, NASTRAN automatically generates the coupled system matrices required by the eigenvalue 
module. 

If the full structure being analyzed exhibits singularities, temporary supports should be specified 
via a SUPORT card. The scalar points that should be supported in this manner are those that represent 
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:t of component rigid-body modes sufficient to define rigid-body motion of the full structure. Note 
t the number of supports used must be equal to the number of rigid-body modes of the full system. 

5 normal choice of temporary support location will be the rigid-body modes of a single component. 

It is recommended that normalization of eigenvectors be done to a unit generalized mass to 
intain a recognizable standard through the transformation back to physical coordinates, 
j Special parameters required for this phase are: 

i NCOMP - The number of components in the synthesis (and for which data have been placed on a 
User Tape). 

PRINT — A positive value will cause the printing of intermediate matrices. 

VECTOR - A positive value will cause the printing of a matrix containing mode vectors in 
component generalized coordinates. 

p partition vectors required are: 

TRUNKi - Identifies modes of component i present on the User Tape that are to be deleted 
! (i.e., not used for the synthesis). M = number of component i modes computed and stored on 
l the User Tape. An element corresponding to a mode to be deleted has value of 1 . 

I POSi - Locates component modes in uncoupled system matrix. M = sum of all component 
! modes to be used in the synthesis. Element has a value of 1 when corresponding to a mode of 
| component i. (The sum of all POSi vectors is a unity vector.) 

LOCi - Locates interface degrees of freedom for component i in truncated uncoupled mode 
vector matrix. M = sum of all interface degrees of freedom. Element has a value of 1 when 
corresponding to a degree of freedom of component i. (The sum of all LOCi vectors is a unity 
vector). 

ferences for MPC, SPC, and eigenvalue METHOD must be present in the Case Control Deck, along 
:h the titles, etc. None of the usual output options is available in this phase. A sample case control 
:k is given in Figure 4. 


RECOVERY OF SYSTEM MODE SHAPES 

e system modal information in physical coordinates can be recovered on a 
nponent-by-component basis using the DMAP sequence of Figure 5, the problem tapes generated in 
ise I, and the User Tape created in Phase II. The User Tape contains not only the modal data but 
> the partitioning vectors supplied in Phase II to define the positioning of the component 
leralized coordinates. These are applied to the system mode vector matrix to truncate it to the 

ieralized coordinates for a selected component, [q^] . The portion of the system mode shape for 

component [</>^pl is expressed in the component global coordinate system by means of 
^-substitution, using the component mode matrix 10^1 > as a transformation: 

r .(i), _ , ,G), , 0), 

[</> ] = [0 ] [q ] (9) 

s c s 

ransient analysis is to be performed, these data are stored on a User Tape for use in Phase IV. Once the 
dal amplitudes of the physical coordinates have been determined, full use of NASTRAN’s output 
ion is available; this option permits printing and/or punching of displacements, constraint forces, 
ernal forces, and internal stresses. 
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REQUEST (UT1,*PF) 

REQUEST! INPT. HO*S.VSN=XXXXX) 

REWIND ( UT1 » INPT ) 

ATTACH(N1S#NASTRAN15» ID=TONE90970) 

MAP (OFF) 

Nib. ATTACH 

CATALOG ( U 1 1 , YOURPERMF ILE » I D=YOURZZ2ZZ ) 

RETURN (INPT) 

exit. 

DMP ( 140000 ) 

RETURN(INPT) 

CATALOG ( UT 1 ♦ YOURPERMF ILE » I D= YOURZZZZZ ) 

(7-8-9) 

ID PHASE . TWO 

APP UISP 

SOL 3.0 

TIME 2 

S( PHASE TWO ALTER PACKAGE GOES HERE) 

CEND 

T l TLE=DEMONSTR AT I ON OF PHASE TWO DECK SET-UP 
MPC = 381 

SPC = 16 

METHOD = S17 
BEGIN BULK 

S (BULK DATA DECK GOES HERE) 

ENDDATA 
( b-7-8-9) 

Figure 4. Typical deck setup for Phase II. 

% phase three dmap SEQUENCE 
begin $ 

INPUTT2 it », »/C»N,-l/V.Y»HOLU=ll 4 
PaRAM //C.N.SUB/V.N.M/V.Y.COMP/ON.l 4 

PaRAM //c , n» mpy/v » n. over/v »N»M/C*N»3 * 

INPUTT2 /CP.RP»TRUNK../V.N.OVER/V.Y.HOLD 4 

INPUTT2 /LAMA'MI.FULL. »/C.N.-5/V»Y»H0LD 4 

PARTN HULL. »CP/. COMPMODE. ./C.N.1/C.N.2/C.N.0/C.N.2 * 

PARTN MODES. TRUNK. /MODEA,. ,/C.N.l/C. N.2/C.N. 2 4 

PARTN QQG. TRUNK »/QA» » */C»N. 1/C . N . 2/C »N» 2 4 

MPYAD MOuEA. COMPMODE. /PHIGG/C.N.O 4 

MPYAD QA, COMPMODE ./QQGG/C.N.O * 

SDR2 CASECC.CSTM.MPT.DIT.EOEXIN.SIL. , .BGPDT.LAMA.OQGG, 
PHIGG. EST . / . OQ , OPHI .0S.0F.PPHI6/C.N.REIG 4 
OKP OPHi »GQ» OF » OS » .//V.N.CARDNO * 

SAVE CAKDuO * 

PARAM //C.N.SUd/V.N.MM/V.N.M/C.N.l 4 
COND OUT 1 , MM * JUMP IF FIRST SUBSTRUCTURE 

INPUTT1 /»»»,/C*N»— 1/V»Y»INP1=1 4 
JUMP GUT2 * 

LABEL OUT 1 * 

OUTPUT 1 . , . , .//C.N.-l/V.Y. INP1 4 

LABEL 0UT2 4 

OUTPUT 1 PHIGG. ., .//V.N.M/V.Y.INPl 4 

END 4 

Figure 5. Phase III DMAP sequence. 
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Since this phase uses a DMAP approach, only output requests have meaning in the Case Control 
Deck and no bulk data are needed. The restart dictionary generated in Phase 1 is used, to which is added 
ja last card reading: 

I n, RE-ENTER AT DMAP SEQUENCE NUMBER 1 

[where n is the appropriate sequence number (see Section 2.2.1 of the NASTRAN User’s Manual). A 
typical deck setup is given in Figure 6. 


REQUEST ( OPTP , HU » S * VSNsYYYYY > 

REQUEST! INPi,HD»S) RING IN 

attach! uti,yourpermfile, id=yourzzzzz> 

REteIND(OPTP,INPl»UTl) 

ATTACH (N15»NASTRANlb»IDsTONE90970) 

MAP (OFF) 

Nl5. ATTACH 
RETURN (OPTP, INP1) 

EXIT. 

DMPU40000) 

RETURN! OP fP, I nPI) 

! 7-8-9) 

ID PHASE, THREE 

APP DMAP 

TIME 

*1 PHASE THREE DMAP SEQUENCE COES HERE) 

S ( PHASE ONE RtSTART DICTIONARY GOES HERE) 

cend 

TiTLE=DEMONGTKATlON OF PHASE THREE DECK SET-UP 
S SPC AND MPC REFERENCES ARE THE SAME AS FOR PHASE ONE 
SPC = 173 
MPC = 4177 
OUTPUT 

SET 16 = 1,3*" THRU 16, 30,35 
SET 31 = l(j6 THRU 135, 167 THRU 194 
VECTOR = ALL 
SPCFORCES = 16 
ELSTKESS = 31 
BEGIN BULK 
EnDDATA 
! (j-7-6-9 ) 


Figure 6. Typical deck setup for Phase III. 

MODAL TRANSIENT RESPONSE ANALYSIS 


iThe transient response analysis via NASTRAN mode synthesis requires user direction to merge the 
required modal data from the component partitions and then perform the response analysis. Since the 
system being analyzed may be quite large, only those degrees of freedom involved in defining external 
loadings or associated with extra points are used. Rigid Format 12, Modal Transient Analysis, 
modified by the alter package of Figure 7, is used. 


The system mode vector matrices for each substructure, written on a User Tape in Phase III, are 
partitioned to extract the portions associated with externally loaded degrees of freedom. If extra 
points are to be added at this stage (e.g., for nonlinear force definitions), the degrees of freedom 
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needed to define the relationship of the extra points to the physical structure must also be included. 
These truncated mode vector sets are merged to form a single truncated system mode matrix. The 
system modes to be used in the response analysis are selected by the user via frequency limits or mode 
count. 

$ PHASE FOUR ALTER PACKAGE 1 

ALTER 23,42 

alter 4a, 91 

$ STORE SUBSTRUCTURE LOAD AND ASSIGN PARTITION VECTORS 
IfoPUTTl /» ,, »/C»N,-l/V»Y,INPl=l * 

0UTPUT1 LOAOlrASSIGNl, , »//V»Y»NCOMP=l/V»Y»INPl * 

0UTPUT2 L0ADi*ASSIGNl»*»//C»N»-l/V»Y»TEMP=i2 * 

$ INSERT ADDITIONAL -OUTPUT 2- STATEMENTS 

* UNTIL all components are treated 

$ temporary disk file (Temp r ia> to reduce tape manipulations 

OUTPUT 1 L0AD2 » ASSIGN2 » » , / /C , N» 0/V » Y » INPl 4 
0UTPUT2 L0AD2# ASSIGN2* , » //C » N» 0/V» Y »TEMP * 
s CONTINUE UNTIL ALL SUBSTRUCTURES HAVE BEEN TREATED 
INPUTT1 /, » * ,/C»N»-l/V» Y> INPl 4 
INPUTT2 /# , , */C»N»-l/V»Y»TEMP 4 
PARAM //C#N*NOP/V»N,TRUE=-l * 

PARAM //C , N » SUB/V » N» K/V » Y » NCOMP/C » N» 1 * 

PURGE PHI A » MAA/TRUE * 

JUMP TOPLOOP 4 

IuPUTTl /PHIA»MAA,GO»GM»KFS/C»N»0/V»Y, INPl 4 DUMMY read 
LABEL TOPLOOP $ 

INPUTT1 /SUBMODE,,, */C»N,0/V,Y, INPl * 

INPUTT2 /LOAD » ASSIGN » » »/C»N»G/V»Y» TEMP 4 

t TRUNCATE MODE VECTOR TO LOADED AND E.P. REFERENCED DOF ONLY. 
PARTN SUbMODE » » LOAD/ » LOADMODE , »/C»N»l/C»N»2/C,N»0/C»N»2 * 

* ASSEMBLE SYSTEM MATRIX 

MERGE, , LOADMODE, ,, ,ASSIGN/LDMD/C,N,i/C,N,2/C,N, 2 4 

AuD PHIA,LDMO/NE*MODE * 

EQUIV NEWMODE»PHIA/TRUE $ 

PARAM //C , N * SUB/ V » N » K/ V » N , K/C » N» 1 4 

COND ENULOOP , K * 

REPT l’OPLOOP ,20 4 

LABEL ENULOOP 4 

CHKPNT PHIA 4 

INPUTT 2 /LAMA , MI , , , /C , N , — 5/V, Y , HOLD— 1 1 4 

CHKPNT LAMA, MI 4 

* MODE MATRIX CONTAINING ONLY NEEDED DOF HAS NOW BEEN FORMED. 

$ CONTINUE STANDARD ANALYSIS 

ALTER 94,94 

alter 98,103 
alter 125 

EQUIV UHVT , UHVTX/NOUE/PPT » PPG/NOUE 4 

COND EOFF , NOUE 4 

PARTN UHVT » * EP/UHVTX » , »/C»N»l/C»N»2/C»N»2 * 

LABEL EOFF 4 

4 OUTPUT SOLUTION SET RESPONSE HISTORY AND LOAD HISTORY 
PARAM //C,N*MPy/V,N,SKIP/V, Y, NCOMP/C, N, 2 4 

OUTPUT 1 UHVTX » PPT , » , //V , N, SKIP/V , Y t INPl 4 

ENDALTER 

Figure 7. Phase IV alter package. 
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Phase IV is similar to Phase II in that only scalar points are used. The scalar points represent all 
aded coordinates and all interface coordinates for extra points, if any. The numbering scheme for 
e scalar points is arbitrary; however, the partitioning vector described below must be compatible with 
e internal sequence of each component, which cannot be altered. Some logical number scheme to 
cihtate bookkeeping is obviously desirable. SEQGP cards may be used to resequence scalar points 
Note that no boundary conditions are specified. They are implicit in the modal information that 
til be used in me analysis. All svstem data are entered in terms of the physical degrees of freedom or 

ore exactly, their scalar point equivalents. NASTRAN automatically performs the required ’ 
msformation to modal coordinates. 

Standard NASTRAN methods for defining external loads sets are used. The loads are applied to 
e scalar points representing the appropriate physical degrees of freedom. It is the user’s 
sponsibihty to ensure that this assignment is correctly specified as no internal check can be made, 
j The interactions between, any extra points and the structure are defined by supplementary 
prices added to structural mass, damping, and stiffness matrices, expanded to accommodate the 
:ditional degrees of freedom. These matrices are input directly via DMIG bulk data cards or formed 
bm transfer functions entered via TF bulk data cards. 

I There is but one special parameter required for this phase: NCOMP - the number of 
bstructures in the system. Partition vectors required include: 

LOADi - Identifies the degrees of freedom in substructure i to be used in the transient 
response analysis. M - g size of the substructure. The element corresponding to the coordinate to 
be used is given a value of 1. Degrees of freedom to be loaded by external forces or that are to be 
referenced by extra points are included. 

ASSIGNi - Identifies the portion of retained degrees of freedom of substructure i in the 
assembled system representation. M = the sum of all coordinates in the structure to be employed 

in the transient analysis. The element corresponding to a coordinate from substructure i has a 
value of 1. 

EP - Identifies position of extra points in augmented modal matrices. M = number of modes 
used in transient analysis plus number of extra points. This vector is required only if extra points 
are used. Further, extra points are always sequenced last in such matrices. A value of 1 indicates 
position of extra point degree of freedom. 

ie Case Control deck includes all references required to perform a transient response analysis. In 
dition, output requests for solution set responses as well as applied and nonlinear force histories are 
ade. A typical deck setup is illustrated in Figure 8. 


RECOVERY OF PHYSICAL RESPONSES 

I 

be response of the entire structural system can be recovered on a substructure-by-substructure basis 
' using the DMAP sequence of Figure 9, the Phase I Problem Tapes, and the Phase IV User Tape. The 
lution set response matrix generated in Phase IV is pre-multiplied by the appropriate partition of the 
stem mode shape matrix to form the substructure response matrix. Then by using standard NASTRAN 
covery methods, the full complement of output options may be exercised. 

This phase is executed in a restart mode; consequently, the mathematical model data are available 
>m the OPTP. (The restart dictionary must be augmented as for Phase III.) The additional 
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REQUEST (INP1 » HD» S» VSNrPPPPP) RJNG in 

ATTACH ( UT1 »YOuRPERMFILE» ID=YOURZZZZZ ) 

RLWIND(INP1»UU) 

ATTACH ( Nib »NASTRANlb» ID=TONE90970 ) 

MAP < OFF) 

Nl5. ATTACH 
RETURN ( I NP1) 

exit. 

DMP ( 140000 ) 

RETURN(INPl) 

(7-8-9) 

10 PHASE# FOUR 
APP U1SP 

SOL 12 » 3 

T X Mf£ 2 

*( PHASE FOUR ALTER PACKAGE GOES HERE) 

TlTLErDEMONSTRATION OF PHASE FOUR DECK SET-UP 
$ K2PP ANO NONLINEAR REFERENCES NOT REQUIRED 
$ FOR LINEAR PRObLEMS. 

K2PP = EXTRA 
NONLINEAR = 5 
DLOAD = 75 
TSTEP = 70 

^ SET 784 = 1 THRU 12»501 THRU 508 
SET 20b = 1 THRU 6 
SET bll - 1 THRU 12 
SET 1800 = 117» 159r 170 » 182» 197 
SET 960 = 960 THRU 969 
NLLQAU = 960 
SDISPLACEMENT = 784 
SVELOC1TY = 208 
S ACCELERATION = 611 
OLOAU = 1800 

begin bulk 

$ (BULK DATA DECK GOES HERE) 

EnDDATA 
(b— 7-8-9) 


Figure 8. Typical deck setup for Phase IV. 

information required consists solely of a partition vector, DISCARD, required to identify the syste 
modes not employed in the transient analysis. M - number of system modes originally oun 
Phase II. Excluded modes are identified by a corresponding element value of 1 . 

Because this phase employs a DMAP approach, no references to boundary conditions, etc., < 
needed. Titling and output requests are chosen as the user wishes for each substructure. Outi 
selections are for physical degrees of freedom and associated functions. A sample deck setup is sho 
in Figure 1 0. 

NASTRAN MODE SYNTHESIS SAMPLE PROBLEM 

The mathematical model of Figure 1 1 has been developed as an aid to demonstrating the procedu 
associated with NASTRAN mode synthesis. While the system is small, it does provide the opportuii 
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* phase five omap sequence 
begin » 

PARAM //C*N*SUB/V»N*MARK/V*Y.COMP=l/C*N*l * 

INPUTTI /**,*/C,N.-l/V*Y.INPl=l % 

INPUTT1 /SUBMODE. * * */V*N»MARK/V, Y» INP1 S 

PARTN SUBMODE* DISCARD* /MODEMaT * » */C*N*l/C*N»2/C*N*2 * 

INPUTTI /UHVTX.PPT* * */C*N*-6/V* Y* INP1 * 

NPYAD MODEMAT»UHVTX*/UDV1T/C,N,0/C*N*1/C*N»1/C,N,1 * 

CASE CASECC * /€ ASExx/C ► N * TRAN/V , N , REPEATr-i/V , N , NOLOOP s 

dpd dynamics»gpl*sil*uset/gpld*sild*usetd»» * * * * * *eqdyn/ 

V * N * LUSET / V * N * LUSETD/ V * N * NOTFL/V * N» NODLT/V * N * NOPSDL/ 

V * N * NOFRL/V * N * NONLFT/V * N * NOTRL/ V * N» NOEED/ 

C * N * 0/ V * N , NOUE * 

VEC USETD/VD/C * N * P/C * N* D/C * N , COMP * 

PARTN UUV1T * *VO/UDSET* * */C*N*l/C*N*2/C*N*2 ‘ 

SDR1 USETD**UDSET***GO*GM**KFS*»/U**QP/C»N*1/C»N»DYNAmICS s 

SDR 2 CASExX*CSTM,MPT,DIT*EQDYN*SILD***BGPDT»PPT,QP,UOVlT,EST,XYCDB/ 

OP * OQ * OU * OS * OF * PUG V/C * N * TR ANRESP * 

CHKPNT PUGV * 

SDR3 OP * OQ • OU * OS * OF * /OP2 * 0Q2 * OU2 * OS2 ► OF2 ♦ ‘ 

OFP 0U2*0P2*0Q2*0F2*0S2*//V*N*CARDN0 $ 

END S 

Figure 9. Phase V DMA P sequence. 

RtGUEST <OPTP*HD,S,VSN=YYYYY> 

RtQUEST( 1NP1*HD*S»VSN=PPPPP) 

Rewind (op tp » inpi ) 

ATTACH (Nib, NASTRAN15*ID=TONE90970) 

MAP (OFF) 

N1 5. ATTACH 

RtTURN(OPTP,lNPl) 

exit. 

DMP ( 140000 ) 

ReTURN(OPTP,InP1) 

(7-8-9) 

ID PHASE* FIVE 
APP UMAP 

TIME 2 

‘(PHASE FIVE UMAP SEQUENCE GOES HERE) 

‘(PHASE ONE RESTART DICTIONARY GOES HERE) 

CtND 

TlTLErDEMONSTKATION OF PHASE FIVE DECK SET-UP 
OUTPUT 

SET 23 = l,3*2b* 16*12 THRU 16 
SET 38 — 38*24,36 
SET 90 = 12*18*24,36*60*66 
SET 91 = 12*18*24*36*60 
SET 92 = 12*18*24*36 
SET 187 = 122*132 
DISPLACEMENT =. 90 
VELOCITY = 91 
ACCELERATION = 92 
SPCFORCES = 187 
ELSTRESS = 23 
ELFOKCE r 38 
BEGIN BULK 

‘(BULK DATA DECK ADDITIONS GO HERE) 

ENDDATA 
( b— 7-8-9 ) 


Figure 1 0. Typical deck setup for Phase V. 
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[dead zone about the null deflection point. Excitation is provided by a sinusoidal force acting vertically 
at the end of the lowest beam. 

For demonstration purposes, the structure is modeled as two components (see Figure 1 1). The 
first, identified as substructure A, consists of the upper and middle beams and the springs between 
i middle beam and the lower beam. The second component, substructure B, is composed of the lower 
beam and its support springs. The component modes for substructure A are developed for the linear 
.portions of the structure only - the nonlinear portions are added when the transient studies are 
performed. 

The first user judgment comes immediately when he must decide which component modes to use 
in the synthesis. For the transient response study to follow, the 30 Hz sinusoidal applied 
: force suggests the desirability of good system definition up to about 1 00 Hz. Based on this objective, a 
I nominal cutoff of 300 Hz for component modes was chosen. For substructure A this presents no 
j problem - the 14 lowest frequency modes plus the three static modes are selected. For 
! substructure B, however, this criterion conflicts with the distribution of modal frequencies. Five 
[vibration modes exist for substructure B that have frequencies under 300 Hz. Adding the three static 
modes gives a total of eight. However, the substructure only has seven independent degrees of 
[freedom. The solution of the conflict was chosen to be a violation of the 300 Hz component mode 
| criterion. For demonstration purposes, still another mode was deleted-, thus the three static modes and 
I three dynamic modes with a highest frequency of 33.7 remain. The resulting system modal frequencies 
i did not show significant error until the 11th elastic mode (70 Hz), where just under 1% error in 
; frequency is observed. However, some degradation of response accuracy results from this 
compromise modeling technique. 

The nonlinear elements between upper and middle beams are integrated with the elastic structure 
definition by addition of extra points and use of a direct input matrix. 

Comparison of the results obtained via mode synthesis with those acquired from the standard 
modal transient study indicates excellent correlation for component A (see Table 1), with errors 
limited to under 1% for displacements and accelerations and under 3% for internal loads. The effect of 
the severe truncation of the substructure B mode set is apparent when comparing results with the 
standard study (Table 1 ). Displacement and acceleration values are quite good, with errors of less than 1%, 
even though the error terms are substantially larger than those for substructure A. Bending moments 
and shears for the selected beam element are in error by 10.7%, because the limited number of elastic 
modes used in the study was inadequate to give good representation. For the more general case, with a 
great many more degrees of freedom in the substructure, it will be possible to employ an adequate 
number of substructure modes to provide sufficient description of elastic behavior. 
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Table 1. Summary of peak responses. 


Response Quantity 

Grid or Ele- 
ment ID 

Quantity 

Standard Analysis 

Synthesized 

Analysis 

Error % 

101 

DISP 

-4.29 

-4.29 

< 0.01 

201 

DISP 

-3.10 

-3.10 

< 0.01 

301 

DISP 

1.11 

1.11 

0.22 

101 

ACC 

5.84 x 10 4 

5.84 x 10* 

0.02 

201 

ACC 

7.12 x 10 4 

7.13 x 10 4 

0.06 

301 

ACC 

1.65 x 10 4 

1.66 x 10 4 

0.5C 

51 

FORCE 

1.43 x 10 4 

1.44 x 10 4 

0.2S 

53 

FORCE 

4.58 x 10 4 

4.58 x 10* 

0 . 0 s 

3 

(bm) a 

1.49 x 10 5 

1.48 x 10 5 

0.2< 

3 

(bm) b 

-6.60 x 10 4 

-6.60 x 10* 

0.0( 

3 

SHEAR 

1.16 x 10 4 

1.15 x 10* 

0.8( 

10 

(bm) a 

4.10 x 10 5 

3.97 x 10 5 

2.9 1 

10 

(BM) b 

4.58 x 10 5 

4.72 x 10° 

2.9 1 

10 

SHEAR 

6.65 x 10 3 

6.78 x 10^ 

1.9 

31 

(BM) a 

0 

0 

£ 


31 

(BM) b 

1.31 x 10 6 

1.17 x 10 b 

A 

10. 

31 

SHEAR 

-5.25 x 10 4 

-4.69 x 10* 

10. 
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